Quantum Magnetomechanics with Levitating Superconducting Microspheres 
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We show that by magneticaUy trapping a superconducting microsphere close to a quantum circuit, 
it is experimentally feasible to perform ground-state cooling and to prepare quantum superpositions 
of the center-of-mass motion of the microsphere. Due to the absence of clamping losses and time 
dependent electromagnetic fields, the mechanical motion of micrometer-sized metallic spheres in the 
Meissner state is predicted to be very well isolated from the environment. Hence, we propose to 
combine the technology of magnetic microtraps and superconducting qubits to bring relatively large 
objects to the quantum regime. 



The field of nano-niechanical resonators aims at cool- 
ing and controlling the mechanical motion of massive ob- 
jects in the quantum regime, and it is a subject of interest 
for both fundamental and applied science The spec- 
tacular progress in this field has led to the achievement of 
ground-state cooling in a high-frequency micromechan- 
ical resonator using cryogenic refrigeration and in 
cavity electro- Q] and optomechanical [1] systems using 
sideband cooling techniques The main challenge in 
these type of experiments is to overcome the heating and 
decoherence produced by clamping losses. These pose an 
even greater challenge for the experimental realization 
of some of the promised applications in the field, such 
as the preparation of quantum superposition states Q, 
which aim at the very fundamental goal of testing the 
validity of quantum mechanics when large masses are in- 
volved Q- Recently, a radical solution to this problem 
has been proposed, namely to unclamp the mechanical 
resonator and to use optical levitation instead [8 12]. 
This is predicted to significantly improve the isolation 
of the mechanical motion from the environment, even at 
room temperature. The archetypical scenario is a dielec- 
tric nanosphere trapped with optical tweezers inside a 
high-Finesse optical cavity. The dependence of the cav- 
ity resonance frequency on the center-of-mass position of 
the sphere yields an optomechanical coupling that can be 
employed to perform ground-state cooling [8|. [p , as well 
as to prepare quantum superposition states (8|, . Re- 
markably, levitation is also the key ingredient in a recent 
proposal to prepare large spatial superpositions, namely 
of the order of the size of the nanosphere 0, [13] ■ 

A common feature in the vast range of optomechan- 
ical systems [ll, including optically levitating mechani- 
cal oscillators |8l4l2|. is that photons are used to cool 
and manipulate the mechanical motion. This introduces 
two important limitations, which are, indeed, the main 
sources of decoherence for optically levitating objects: 
(i) scattering of photons, that produces position localiza- 
tion decoherence @, i. Ill, and ( ii) absorption of pho- 
tons, that increases the bulk temperature of the object, 
and, consequently, also the decoherence due to emission 
of black body radiation 0-0, In this Letter, we 



propose a levitating mechanical resonator experimental 
setup which does not employ photons but magnetostatic 
fields. Therefore, it does not only avoid clamping losses 
due to levitation, but it also circumvents limitations (i) 
and (ii) due to the absence of photons when trapping, 
cooling, and coherently manipulating the object. Hence, 
this mechanical oscillator is predicted to be very well iso- 
lated from the environment. This leads to very large me- 
chanical quality factors and long coherences times (due 
to negligible black-body radiation), the latter being spe- 
cially relevant in time-of- flight experiments 0, El . The 
proposal consists in magnetically trapping a supercon- 
ducting microsphere in the proximity of a quantum cir- 
cuit [e.g. LC oscillator or fiux qubit) [l4|. The mag- 
netic field employed in the magnetic trap is expelled (see 
Fig. [1]), due to the Meissner effect, from the supercon- 
ducting sphere. The flux passing through the pick-up coil 
of the quantum circuit depends on the center-of-mass po- 
sition of the sphere. This leads to a significant quantum 
magnetomechanical coupling with the center-of-mass mo- 
tion. Here we show that ground-state cooling and the 
preparation of quantum superpositions can be achieved, 
within the same experimental setup and with feasible pa- 
rameters, for superconducting spheres {e.g. Pb) in the 
/im regime with masses of ^ 10^'* amu. 

More specifically, we consider a superconducting mi- 
crosphere of radius R and mass M that is in the Meissner 
state, cooled below a certain temperature Tq. The pen- 
etration length A and coherence distance ^ are such that 
i? 3> A,^. In this regime, one can approximate that the 
magnetic induction B is zero in the whole interior of the 
superconductor. The microsphere is confined into a 3D 
harmonic potential created by a magnetic microtrap 15]. 
While the proposal does not rely on the particular trap- 
ping scheme, we consider here a quadrupole trap created 
by two circular coils of radius ^ 3> i? in an anti-Helmholtz 
configuration (AHC), i.e. the coils are coaxials, separated 
by a distance Z, and with opposite current intensity /, 
see Fig[TJ This trap creates an harmonic potential of the 
form T4rap = M[ujIx^ + uj'j_{y^ + z^)]/2. The trapping 
frequency cot can be obtained analytically using the im- 
age method due to the cylindrical symmetry of the 
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FIG. 1: (color online) (a) Scheme of the AHC with a pick- 
up coil, (b) Quadrupole field created by the AHC. (c) Field 
expelled by the sphere in the Meissner state. 



system, see the Supplemental Material (SM). The trans- 
verse frequency can be derived by replacing the sphere 
by an effective magnetic moment. Their expressions are 
given by wt — 1.05-\/(/io/p)///^ and uj±_ = LUt/2. fiQ is 
the vacuum permeability and p the density of the micro- 
sphere which is assumed to be homogeneous. An impor- 
tant remark is that the field at any point of the sphere is 
required to be smaller than the critical field Bq in order 
to allow superconductivity; this yields an upper bound on 
the radius of the sphere R < -Rmax — 0.98i?c / (wt y/f^op) ■ 

At some distance d from the center of the trap and 
coaxially with the AHC (Fig [1]) a pick-up coil of ra- 
dius r is placed (the perpendicular configuration might 
also be considered, see SM). The pick-up coil is c on- 



nected to a quantum LC oscillator or to a fiux qubit 14 1 
placed outside the magnetic field of the trap. As shown 
below, the quantum magnetomechanical coupling for 
both cases depends on the dimensionless parameter rj = 
Xzp<i>^xj(0)/<I>o, where x^p — \/h/ (2Mwt) is the mechan- 
ical zero point motion, $o = T^h/e is the flux quan- 
tum, and $J,xt(0) is the derivative with respect to the 
axial motion of the center-of-mass evaluated at its equi- 
librium position of the fiux threading the pick-up coil. 
^cxt (0) can also be evaluated analytically using the image 
method [l6j and its approximated expression for R <^l,r 
is given by ^'^^^{Q) w 2.7 /l'^)R^r^ / {<P + r^f^ (see 
SM for further details). 

For an LC oscillator, consisting of an inductor L con- 
nected to a capacitor C, the Hamiltonian is given by 
-ffLC = [* - $cxt(i)] V(2i) + QV(2C), with [$, Q\ = ih. 
By expanding $0x1(2;) linearly in x, one obtains the lin- 
ear Hamiltonian Hi^c — fiujhco-'' a + hgi,c{a^ + o){o^ + 
We used x = Xzp{i>^ + b) and $ ~ $zp(a^ + a), where 
yjh/ {2CijJ-lg) is the zero point flux, wlc = 1 /VLC 



is the LC resonator frequency, and ^lc — ^hcVi with 
ELC = ^o^zp/ {hL), is the magnetomechanical coupling 
rate. Alternati vely , one can consider a three junction 
flux qubit [13, [13 ■ The two qubit states correspond 



zp 



to persistent currents of amplitude Ip flowing clockwise 
or counterclockwise in a loop of radius r. When the 
qubit is operated in the vicinity of the degeneracy point 
/(<I>oxt) = <J'cxt/*I>o - 1/2 « 0, the Hamiltonian of the 
qubit in the basis of the persistent current states reads 
Hg = —heaz/2 — hAcTx/2. Here (t^ {i — x,y,z) are the 
usual Pauli matrices, e — J^/($oxt)/^ is the bias (where 
v = 2'I>o^p), and A is the tunneling amplitude. Thus, by 
expanding e{x) ~ e(0) -|~e'(0)x to flrst order in x, one ob- 
tains the quantum magnetomechanical Hamiltonian (the 
harmonic oscillator energy is also included). 



H 



MM 



^ = ujtb^'b - ^a;, - - gad-z{b^ +b). (1) 



Here, e = e(0), and go 
tomechanical coupling. 



~ vrj is the non-linear magne- 
The energy levels of the qubit 



are separated by a frequency Wg = v -I- A^ . It is re- 
markable that with typical numbers, see below, v is only 
two or three orders of magnetic lower than clc since the 
coupling to the flux qubit provides a non-linearity to the 
mechanical motion. For this reason, we concentrate here- 
after on the coupling with the fiux qubit in order to dis- 
cuss ground-state cooling and the preparation of quan- 
tum superpositions. However, due to its larger coupling, 
the LC oscillator would be definitely more appropriate 
for performing cooling or any other Gaussian dynamics 
{e.g. preparation of squeezed states or entanglement and 
teleportation in the case of several spheres) . In this case 
a parametric driving in the magnetomechanical coupling 
{e.g. with a time-dependent inductance) would be re- 
quired. We remark that while we do not consider it here, 
coupling to the transversal motion also leads to a signif- 
icant magnetomechanical coupling. 

Decoherence in superconducting qubits is modeled 
with a master equation containing the following Lind- 
bland terms [3] (written in the eigenbasis of the qubit): 
spontaneous emission £o[p] = ro(2o'_/5(7+ — pa^a^ — 
d+(j^p)/2, and pure dephasing F^, L^[p\ = T^{dzP^z - 
p)/2. The decoherence rates are related to the relaxation 
time Ti = I/Fq, which is the time required for the qubit 
to relax from the excited state to the ground state, and to 
the dephasing time T2 = (Fo/2-f F^)"^, which is the time 
over which the phase difference between two eigenstates 
becomes randomized. 

The distinctive feature of this proposal is that the de- 
coherence in the mechanical oscillator, when accounting 
for known sources, is predicted to be very small since 
both (i) clamping losses and (ii) scattering of photons 
are absent. Similarly as in the proposal of optically levi- 
tating dielectrics Ell ' other sources are also negligi- 
ble: (iii) damping created by the background gas yields a 
mechanical quality factor that can reach extremely high 
values at sufficiently low pressure, this is given by 

Qair = t^t/7air > lO", where 7air = l&P/{T:vRp), 

P ~ 10~^° Torr is the environmental pressure, and v 
the thermal mean velocity of the air molecules; (iv) de- 
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coherence due to blackbody radiation 

0, S m yields 

rates at the Hz regime due to the cryogenic bulk temper- 
atures of the metallic superconducting microsphere (in 
optical levitation the bulk temperature of the object is 
heated due to light absorption); (v) internal vibrational 
modes are decoupled (due to their higher frequency) to 
the center-of-mass motion for micrometer-sized objects, 
see [ill for a detailed analysis based on quantum elastic- 
ity. Other sources of decoherence particular to this pro- 
posal are: (vi) damping due to hystereses losses in the 
superconducting coils yield Q ^ 10^° (specially for small 
fluctuations in the position, see SM for details) ; (vi) fluc- 
tuations in the trap frequency leads to decoherence with 
a rate given by = ttw^ S'i^(2a;t)/16, where Sui{'iiOt) 

is the one-sided power spectrum of the fractional fluc- 
tuation in the resonance frequency (see SM for details). 

~ Hz can be obtained for yj S^{2iiJt) = 10^^/\/Hz; 
(vii) fluctuations in the trap center also lead to deco- 
herence with the rate Tx = ttw^ S'x(wt)/4a;2p, where 
in this case Sx is the one-sided power spectrum of the 
position fluctuations, ^ Hz can be obtained for po- 
sition stability y/S^{uJt) / Xzp ~ 10~'*/VHz; (viii) within 
trapping frequencies in the kHz regime, which are much 
smaller than the energy gap (100 GHz), the supercon- 
ductor can be considered to act instantaneously to ex- 
ternal fields; (ix) and finally, we remark that the flux 
qubit present in the setup can always be decoupled from 
the center-of-mass motion by switching off the driving 
such the coupling is off-resonant (see below). Thus, at 
sufficiently low pressure and for superconductors in the 
Meissner state, the center-of-mass of micrometer-sized 
metallic spheres is (according to the sources of decoher- 
ence that we have considered) effectively decoupled from 
the environment. This assumes very stable traps (both 
in the frequency and in the equilibrium position). In the 
actual experiment, other uncontrolled sources of deco- 
herence might be relevant, for instance, coupling to bond 
paramagnetic centers on silicon surface have been exper- 
imentally observed at distances of few micrometers [20|. 

We now focus on ground-state cooling of the mechani- 
cal motion. To this end, a resonant coupling between the 
qubit and the mechanical resonator is required [2ll . 
Note however that the energy splitting of the qubit Wg (in 
the GHz regime) is much larger the mechanical frequency 
ujt (in the kHz-MHz regime) of the oscillator. This sce- 
nario has been studied both experimentally and theoret- 
ically in 18| , where a flux qubit has been coherently cou- 
pled to a slow LC oscillator. This is achieved by driving 
the flux qubit with an applied ac flux with frequency Wd 
(in the GHz regime) and amplitude il (in the kHz-MHz 
regime). In this case, the total Hamiltonian reads Ht = 
HuM+Hdi-ivc, where i/drive — hncoa{ujdt)az- Recall that 
the dynamics of the qubit and the mechanical oscillator 
(which is decoupled from the environment) is given by the 
master equation p = —i[H^,p\/h + LQ[p\ + L^[p\, where 

is the total Hamiltonian written in the eigenbasis of 
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FIG. 2: (color online) (a) The steady state phonon number 
occupation (n)ss (assuming Text = 0) and (b) the cooling rate 
r over r = cos^(Q:)/ro are plotted as a function of j3 for 
r^/Fo = (solid blue line), V^/Vq = 0.1 (dashed red line), 
and r^/Fo = 1 (dotted orange line). 



the qubit. By moving to a rotating frame of the qubit at 
frequency Wd, performing a rotating wave approximation 
(RWA) (valid provided that Wd ^ i^s ^ "^tjffo), trans- 
forming to the diagonal basis of the qubit, moving to the 
interaction picture, and performing a second RWA (as- 
suming ujt ~ Cbs^ go, Tq, Fip, where = {Suj"^ + fl'^Y^''^, 
Suj = Wd— Wg, Q = Qsina, and tana ~ A/e), one arrives 
at 



g(CT_6t + H.c.),p +Cr[p 



(2) 



The effective magnetomechanical coupling is given by 
g = 5o cos a sin /3, where tan/3 — fl/Suj. The dissipa- 
tion of the qubit is given by Ct[p] and contains both 
dephasing r'^{azP<7z — /5)/2, and transitions in both di- 
rections Fj^(-|-) (2(Tzpp(7± — pa±a^ — (t±(t^p) /2. The rates 
are given by F* ~ T^p cos^ /? + Fq sin^(/3)/2, and F^(.|.) — 
sin2(/3) -I- Fo(l ± cos/3)V2. Typically ~g < T^^, F^, 
which allows to adiabatically eliminate the qubit |23i] . 
This leads to an effective master equation describing the 
mechanical oscillator density matrix that can be used to 
obtain a dynamical equation for the mean phonon num- 
ber occupation fi = wh, namely {n) = — F(n) +A-|_-|-rcxt. 
The cooling rate is defined as F = — where 



A+ = 



2g' 



(2F*-fFt + F;) 



It 



(3) 



The rate Foxt takes into account external sources of heat- 
ing in the mechanical oscillator {e.g. due to trap fluc- 
tuations) which are assumed to be much smaller than 
F. The final phonon number occupation is given by 
(7i)ss = (A+ +rcxt)/r. In Fig.[2]it is shown that ground- 
state cooling can be achieved within a wide range of 
(3 with cooling rates of the order of ~ 5^ cos^(a)/ro 
(a < 7r/2 is also required). 

The non-linearity given by the presence of the qubit 
can be readily employed to prepare quantum super- 



positions of the mechanical oscillator [2l|. First, we 
rotate the Hamiltonian Eq. ([1]) to the eigenbasis of 
the qubit, we move to the interaction picture (using 
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FIG. 3: (color on-line) Scheme of the protocol to prepare 
quantum superpositions of a mechanical oscillator using the 
parametric coupling to a qubit. (a) The state |-l-)|0) that is 
prepared t — and recovered ai t = 2t* . (b) The superposi- 
tion state Eq. Q that is created at t — t* . 



as a free part the qubit energy term), and perform 
the RWA (valid provided cjg ^ go)- We arrive at 
^iiM = + fi^sCTz - g^zx/x^p, where g = go cos a 
and = p2/(2M) -f This can be re-written 

as -ff^M = ^s/2o-z + T\x^z)H^f{xcrz), where r(a) = 
exp[— ipaxzp/ft] is the translational operator (such that 
T'^ {a)xT{a) — x + ax^p), X = "^g/^t is a dimension- 
less parameter, and the constant term — fig^/cjtl has 
been dropped out. The Hamiltonian written in this 
form points out the key property used in the protocol, 
that is, that the center of the harmonic trap depends 
on the spin state of the qubit, see Fig. EJa). Hence, 
by initially preparing the joint system into the state 
l*o) = l+,0) = (I t,0) + I ;,0))/%/2, that is, the qubit 
in the superposition state |+) = (| t) + I i))/V^ and the 
mechanical oscillator is in a pure state |0) with (x^) — a^, 
the joint state evolves after time — 7r/wt into the en- 
tangled state 



1 

V2 



T(-2x)|t,0)+T(2x)U,0) 



(4) 



Hence, the mechanical oscillator is in a spatial quan- 
tum superposition state (Fig. [2Ib)) where the wave 
packets are separated by a distance Ig — ^x^zp = 
Sx^pg/uf The overlap is given by (0|r^(-2x)T(2x)|0) = 
exp [— /(8ct^)] , and thus 8cr^ < ll is required. Due to 
levitation, this challenging condition [2l| can be guaran- 
teed by initially squeezing the ground state of the me- 
chanical oscillator, see SM. The superposition can be 
probed by performing tomography of the qubit during 
the evolution of the joint state within a time window 
t € [0, 2t*], such that the collapse at t = t'^(due to the 
entangled state Eq. (|4])) and revival at i = 2t* (due to 
the product state) of the purity of the qubit state can 
be observed. Decoherence in the qubit can be neglected 
provided 2t* <C T2. Note that for long coherence times 
in the qubit, the superposition size could be increased by 
performing a spin flip after each ti, evolution, or alterna- 
tively, by opening the trap to a frequency w' <C wt- 

Regarding the experimental feasibility of this proposal, 
in the following we choose some set of parameters that 



fulfill the conditions that are required. We consider a 
microsphere made of Pb, with density p = 11360 Kg/m'^, 
penetration depth A = 30.5 nm and coherence length 
^ = 96 nm (evaluated at T = 0), critical temperature 
Tc = 7.2 K, and critical field Bq = 0.08 T. The mag- 
netic trap is made with AH coils of radius I = 25 fj,m and 
current / = 10 A (assuming materials with critical cur- 
rent Jc = 7 X 10"A/m2 gj]) such that ut ~ 2tt x 28 kHz 
and i?max = 3.7 /im. Thus, we consider a sphere of radius 
R — 2 /zm. The pick-up coil has a radius r = 24.5 /im 
and is placed at d = 17.5 /zm (outside the AHC). This 
yields rj = 1.3 x 10~^. For the flux qubit we use [l3| 
= A = 27r X 10 GHz and Ti = T2 = 10 ^is, and 
hence we obtain g^ — 2tt x 1.3 kHz, Fq = 27r x 16 kHz 
and F^ = Fo/2. For an LC oscillator with C = 1 



pF and L = 0.1 nH [lj|, the magnetomechanical cou- 



pling is nearly two orders of magnitude larger, namely 
.9lc = 27r x 93 kHz. This set of parameters allows for 
ground-state cooling and the preparation of quantum su- 
perposition states using the squeezing of the ground-state 
wavefunction. 

We have shown that micrometer-sized superconducting 
metallic spheres containing ~ 10^^ atoms can be cooled 
down to the ground state and prepared into superposi- 
tion states. From a broader perspective, while cavity op- 
tomechanics is using the technology developed for laser 
trapping, cooling, and manipulation of atoms and ions, 
here we propose to merge the technology of the growing 
fields of magnetic trapping of atoms and superconducting 
qubits to bring massive objects into the quantum regime. 
The combined properties of levitation, low bulk temper- 
atures, and large masses, makes this setup ideally suited 
to design and implement protocols where the objects is 
released from the trap in order to expand the wavefunc- 
tion. This can be used to create large superpositions fl^l 
in order to test fundamental questions ^7|] and to design 
ultra-high sensitive detectors. 
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